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Abstract 

We compute the total cross-section for Higgs boson production in bottom-quark fusion using 
the so-called FONLL method for the matching of a scheme in which the 6-quark is treated as 
a massless parton to that in which it is treated as a massive hnal-state particle. We discuss 
the general framework for the application of the FONLL method to this process, and then 
we present explicit expressions for the case in which the next-to-next-to-leading-log hve-flavor 
scheme result is combined with the leading-order 0{a‘^) four-flavor scheme computation. We 
compare our results in this case to the four-and five-flavor scheme computations, and to the 
so-called Santander matching. 
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In perturbative QCD, processes involving bottom quarks can be computed within different 
factorization schemes. One possibility is to use a five-flavor, or massless, scheme, in which the 
6-quark is treated as a massless parton. In this scheme, collinear logarithms of fijp/ml (with 
fxp the factorization scale) are resummed through QCD evolution equations, but corrections 
suppressed by powers of ml /are neglected. Alternatively, one may use a four-flavor, massive, 
or decoupling scheme, in which the 6-quark is treated as a massive particle, which decouples from 
evolution equations and the running of a*, but full dependence on mb is retained. Generally, 
of course, results in the two scheme may differ by a large amount: indeed, the leading-order 
predictions for Higgs boson in bottom-quark fusion HHa may differ by up to one order of 
magnitude [5], though the disagreement is reduced if the factorization and renormalization 
scales are chosen to be smaller than mn (which may well [GllTO] be more appropriate) and 
higher perturbative orders are included. 

The five-flavor scheme is more accurate for scales /i^ ^ m^, while the four-flavor scheme is 
more accurate close to threshold, though of course if the four-flavor computation is performed 
to high enough order in perturbation theory it will reproduce the five-flavor scheme result (the 
converse is not true, because mass corrections are not included in the five-flavor scheme at any 
perturbative order). It is therefore advantageous to combine the two computations into one 
which is accurate at all scales. A phenomenological way of doing so, the so-called Santander 
matching, has been proposed in Ref. it consists of simply interpolating between the four- 
and five-flavor scheme results by mean of a weighted average, such that in the two limits 
p/mfe ^ 1 or /i/mb ~ 1 the massless or massive results are respectively reproduced. 

However, a more systematic approach which preserves the perturbative accuracy of both 
computations may be desirable. One such approach, the FONLL method, was proposed in 
Ref. [12] in the context of hadro-production of heavy quarks, and extended to deep-inelastic 
scattering in Ref. [T3|. The basic idea of this method is to expand out the five-fiavor-scheme 
computation in powers of the strong coupling as, and replace a finite number of terms with their 
massive-scheme counterparts. The result then retains the accuracy of both ingredients: at the 
massive level, the fixed-order accuracy corresponding to the number of massive orders which 
have been included (FO, or fixed order), and at the massless level, the logarithmic accuracy of 
the starting five-flavor scheme computation (NLL, or generally subleading logarithmicQ) . 

It is the purpose of this paper to present the application of the FONLL scheme to Higgs 
production in bottom-quark fusion, focusing for definiteness on the total cross-section. In the 
rest of this paper we will follow the notation and conventions of Ref. [13] . 

The total cross-section a in the five-flavor scheme has the form 

dxidx2'^ fP {xi, alf {xi,X2,af\fx^)) , ( 1 ) 

ij 

where the sum runs over the 10 quarks and antiquarks and the gluon, and the 6 quark and 
antiquark are treated as the other partons, which in particular contribute to the running of 
For simplicity we omit the dependence of the hard cross-section on the renormalization 

^We will consistently use the notation N^LL to refer to the resummation of collinear logs of the heavy quark 
mass, i.e. by LL we mean a computation in which Tas In is treated as order one (0;°). 
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and factorization scales, which henceforth we will assnme to be chosen eqnal to hr = = fji-, 

unless otherwise stated. 

In the fonr-flavor scheme it has the form 




JJ dxidx2 f, 

ij 


( 4 ) 




a. 


( 4 ) 


(4) 

Xi,X2, —,al ’ 


mt 



( 2 ) 


where now the sum only runs over the four lightest quarks and antiquarks and the gluon, the 
6-quark decouples from the running of and the DGLAP evolution equations satisfied by 
{xi, but full rrih dependence of the partonic cross-section is retained. 

In order to carry out the FONLL procedure, we need to express the four-flavor scheme 
cross-section, Eq. ([2]), in terms of and so that their perturbative expansions can be 
compared directly. The coupling constant and the PDFs are related in the two schemes by 
equations of the form 


i=2 


(3) 


jX y ^ \y / 

where 

L = In (5) 

and the sum runs over the eight lightest flavors, antiflavors, and the gluon, while the index 
i takes value over all ten quarks and antiquarks and the gluon. The coefficients Cj(L) are 
polynomials in L, and the functions Kij can be expressed as an expansion in powers of with 
coefficients that are polynomials in L. 

The hrst nine equations (jl|) relate the eight lightest quarks and the gluon in the two schemes 
and can be inverted to express the four-flavor-scheme PDFs in terms of the hve-flavor-scheme 
ones. The last two equations, assuming that the bottom quark is generated by radiation from 
the gluon (i.e. no “intrinsic” HI bottom component) express the bottom and anti-bottom 
PDFs in terms of the other ones. In particular, this assumption implies that the b quark and 
antiquark PDFs are equal to each other, fjf"' = Inverting Eqs. (jSHl]) and substituting in 

Eq. ([2]) one can obtain an expression of in terms of and 




v=q,9 


dxidx2 ^ jf\xi,fi‘^)ff\x2,n^)B^f lxi,X2,-^,af\fi^) 

\ TTlu 


( 6 ) 


where the coefficient functions Bij are such that substituting the matching relations Eqs. fl^ -fHl) 
in Eq. ([HD the original expression Eq. ([2]) is recovered. Note that in the course of the procedure 
of expressing in terms of and , subleading terms are introduced, because Eqs. ([3}11]) 
are only inverted to hnite perturbative accuracy. It follows that the expressions Eq. ([2]) and 
Eq. ([6]) of actually differ by subleading terms. Henceforth, for we will use the expression 
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Eq. ([6]), and avoid any further reference to and therefore, from now on as and fi will 
denote the hve-flavor scheme expressions. 

In order to match the two expressions for a in the hve-flavor scheme, Eq. o. and in the 
four-havor scheme, Eq. (E]), we now work out their perturbative expansion. Using DGLAP 
evolution, the 6-PDF, can be determined in terms of the gluon and the light-quark 

parton distributions ff'^ at the scale convoluted with coefficient functions expressed as 
a power series in ai^\ with coefficients that are polynomials in L. The hve-havor-scheme 
expression Eq. ([T]) may thus be written entirely in terms of light-quark and gluon PDFs: 


dxidx2 ^ f^^\xi,n^)ff{x2,fJ,^)A[f(xi,X2,L,af\fi^)), (7) 

f5) 

where the Aij coefficient functions are given by a perturbative expansion of the form 

N OO 

A^f {xi,X 2 , L,af\fj.^)) = {af\i?)Lf , (8) 

p=0 k=0 

with at leading order iV = 0, and at N™LO order N = m. 

On the other hand, the four-havor-scheme expression Eq. ([6)), as mentioned, is also written 
in terms of the light-quark PDFs, with coefficient functions Bij which can also be expanded in 

t ( 5 ) 
power of Qfs , 




(5) 

Xi,X2, > 

mi 



N 





(9) 


where N is the order of the expansion needed to reach the desired accuracy. It follows that the 
sum of all contributions to the four-havor-scheme expression Eq. ([9]) which do not vanish when 
jj? S> ml must also be present in the hve-havor-scheme result. 

These contributions provide the massless limit of in the sense that 


lim 

>-0 




o(0),(p) 

ij 



0 . 


( 10 ) 


In other words, is obtained from B^^ by retaining all logarithms and constant terms 

and dropping all terms suppressed by powers of m^//i. Given that these terms are also present 
in the hve-havor-scheme calculation, we can also write 


and 


B, 


(o).(p) 




( 

Xi,X2, 

V 



{p-k),(k) 

ij 


{Xi,X2)L^ 


B. 


( 0 ) 


1-3 


Xi,X2, ’ 

ml 


N 


( 9 ‘^))=E( 

p=0 


a 


( 5 ) 


(/^^)) 


P r(0),(p) 
ij 


, X2 1 2 

mi 


(11) 


( 12 ) 
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We finally define the massless limit of the fonr-fiavor-scheme cross-section, namely 


a 


(4).(0) _ 


JJ dxidx2 ^ /, 




( 0 ) 

ij 




W (5) 

Xi,X2, —,al 
mf 


(p-) 


( 13 ) 


The FONLL method can thus be stated as follows: replace in the five-flavor scheme expres¬ 
sion, Eq. d^, all contributions to the expansion Eq. (|8]) of the coefficients (^xi,X 2 , L, 

which appear in Eq. flTT]) . with their fully massive expression from Eq. ([9]). In this 

way, all mass suppressed effects that are not present in Eq. ([T]) but are known from Eq. (|2]), 
are included. More symbolically 

^FOVLL^^(4)^^(5)_^(4),(0)_ (^4) 

If the five-flavor scheme computation is performed to N^LL accuracy, and the replacement is 
performed up to fixed N-’LO in aW, the final result retains N'^LL accuracy at the massless level, 
and N-^LO accuracy at the massive level. 

In Ref. [13], three combinations were considered specifically in the case of deep-inelastic 
scattering: namely FONLL-A, corresponding to NLL-LO, FONLL-B, NLL-NLO, and FONLL- 
C, NNLL-NLO (where by “leading” we always mean the first order at which the result does 
not vanish, assuming no intrinsic heavy quarks). In deep-inelastic scattering, the leading order 
is 0{a^) (parton model) in the five-flavor scheme, and 0{as) in the four-flavor scheme: there 
is thus a mismatch by one order, and therefore FONLL-A is the simplest nontrivial scheme. In 
the case of Higgs production in bottom fusion, the mismatch is now by two orders: the leading 
order is 0(a°) (parton model) in the five-flavor scheme, and 0{a‘^) in the four-flavor scheme. 
The simplest nontrivial case, which we will also refer to as FONLL-A, is thus NNLL-LO; we will 
then call FONLL-B the NNLL-NLO combination and FONLL-C N^LL-NLO. In the five-flavor 
scheme, the result is known up to NNLO na, thereby allowing for an NNLL computation when 
used in conjunction with NNLO PDFs, and in the four-flavor scheme up to NLO [iGlfTT] . hence 
in principle FONLL-A and FONLL-B are accessible using current knowledge. 

We now work out Em (fT4|) explicitly for Higgs production in bottom-quark fusion, in the 
simplest FONLL-A caseo To NNLL, the partonic cross-section must be computed up to order 
(P(a^): it then receives contributions from the following sub-processes: 

• 0{l)^bb^h 

• 0{as) bb ^ h (1-loop), bg —)■ hb, bb —)■ hg 

• 0{al) bb ^ h (2-loop), bg —)■ hb (1-loop), bb hg (1-loop), bq —)■ hbq, gg hbb, 
bb hbb, qq hbb. 

matched computation for the related process of Higgs production in top fusion has been presented 
recently [TS], based on a modified version of the ACOT [12] matching scheme, which for NLO deep-inelastic 
scattering is known [20] to coincide with FONLL-A; however, in this work only terms up to NLL in the five-flavor 
computations are included. 
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Figure 1; Leading-order (a) and next-to-leading order (b-c) contributions to the hard cross- 
section in the hve-flavor scheme. To order 0{al) these processes receive 2-loop corrections (a) 
and 1-loop corrections (b) and (c), respectively. 


The LO diagrams are shown in Fig.[T] The full calculation up to can be found in Ref. [15] . 

The relevant perturbative orders in each parton channel are thus 

{xi, X2, = a^'^’^^\xi,X2) + af\fi^)a^'’’^^\xi,X2) 

+ + 0{al), (15) 


^(5) 

(a^i, 

X2, 



.X 2 ) 










[xi, 

X 2 ) 

+ 


(16) 

^(5) 

{xi, 

X2, 

«P(y)) 


{xi, 

X 2 ) 

+ 


( 17 ) 

<y(5) 

99 

{xi, 

X2, 



(a^i, 

X 2 ) 

+ 


(18) 

^(5) 

{xi, 

X2, 



{xi, 

X 2 ) 

+ 


(19) 

-(5) 

^qq 

{xi, 

X2, 



{xi, 

X 2 ) 

+ 

0{al). 

(20) 


In the four-flavor scheme, the LO 0(ag) result corresponds to the gg —> hbb and qq hbh 
sub-processes shown in Fig. O The computation of this process in the four-flavor scheme is 
formally identical to that of associate production of a Higgs boson with a ti pair, hrst performed 

in Ref. [1]. 

We can now match the two expressions. First, we note that in the FONLL-A scheme the 
four-flavor scheme result is included to lowest nontrivial order: therefore, we can simply replace 
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Figure 2: Leading-order order contributions to the four-flavor scheme. Not shown are diagrams 
that can be obtained by crossing the initial state gluons, or radiating the Higgs off an anti 
bottom quark. 


in it and with their five-flavor scheme counterparts, as the difference is higher order in 
as and thus subleading. We thus simply have 

Bij X2, = d--f £|, Q!s(/i^)^ 0{al). (21) 

We also need the massless limit of the four-flavor scheme result: recalling that it starts at 
order a^, and using the general expressions Eqs. flTTll - ffT^ , we conclude that it must have the 
form 

Blf{xi,X2,L,as = {asf {xi,X2,L) + O(a^) 

= {asf (^Alf’^^\xi,X2) + A[f'^^\xi,X2)L + Alf’^‘^\xi,X2)L‘^^ + 0{al). 

(22) 

The easiest way of determining the coefficients is to start with the five-flavor scheme 

expression Eq. ([1]) and expand the bottom PDF in power of as, 

d A A + 0(aj). (23) 

V X y \ y / 

where 

P„(v) = Tb. [y + (1 - !/)T . (24) 
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We get 


Af^’^^\xi,X2) = af^'^‘^\xi,X2), 

Afl^^\xi,X2) = dyPqgiy) (dlf^^\xi,yx2)+ alf^^\yxi,X2)^ + (&->&), 

Af^’^^\xi,X2) = -7^-^ [ [ dyidy2Pqg{yi)Pqg{y2)dlf’^^\yiXi,y2X2) + {b ^ b), 

l^TTj Jq Jq 

SO that 


(xi,X2,L,a,) = (11,12) + A‘]’’’‘’(li,l 2 )i + 


(25) 

(26) 

(27) 

(28) 

(29) 


We now have all the ingredients which enter the FONLL-A expression, 
purposes, we introduce a formal expansion of the cross-section of the form 

^FONLL-A ^ ^FONLL-A,(0) ^ ^^2yFONLL-A,(l) ^ 2 ^fONLL-A,(2) 


For book-keeping 
+ 0(af). (30) 


where it is understood that only the coefficient functions b\^\ aIj'^ and in Eqs. (jb]), ([7]) 
and (1121) respectively are expanded, but not the PDFs. The expansion is formal in that, as we 
have just seen, the nominally 0(a°) contribution really starts at O^a"^) once one substitutes 
the explicit expression Eq. fl23l) of the 6-quark distribution, as it should be in order for it to 
match the four-flavor scheme expression. 

Be that as it may, since the four-flavor scheme starts at 0(ag), the hrst 

two terms in the expansion Eq. fl30l) coincide with the hve-flavor scheme expressions: 


^FONLL-A,(0) 


dxidx2fb{xi, ix^)f-b{x2, {xi,X2) 


(31) 


0-FONLL A,(1) ^ JJ dxidx2^fb{xi,i/)fi{x2,y‘^)dl^’^^^{xi,X2) 

+dfb’^^\xi,X2) [{fg{xi,y'^)fb{x2,y^) + {xi X 2 )) + (6 -A 6)] |. (32) 

The 0(ag) contribution can be written as the sum of two terms: four-flavor scheme, and 
difference between the hve-flavor and the massless limit of the four-havor scheme. The former 
is simply given by the leading-order partonic cross-section in the four-havor scheme.The latter 
is given by 

^FONLL-A,(2)^^(4),(2)^^(d),(2)^ ^ 33 ) 

where 


^{d),{2) ^ ^(5),(2) _ ^(4),(0),(2) 


(34) 


a 


(4),(0),(2) _ 


dxidx2 ^ {xi,X2,L,as). 


(35) 


ij=q,9 


and 



We get 


^{d),(2) = / / dXidX2 < fb{Xi,fi‘^)fl{x2,fi'^)dl^’^‘^\xi,X2) + 


+ fb{xi, fX^)fb{x2, ld^)d^^’^‘^\xi, X2) + 

+ afJ’^^\xi,X 2 ) [{fg{xi, ^‘^)fb{x2, fi^) + {xi ^ X2)) + {b ^b)] + 


a 


( 5 ),( 2 ) 
qb 


(a^i, X2) [{fg{xi,fi'^)fb{x2, /i^) + {xi X2)) + {b b,q ^ q)] + 


2tt 


dyPqgiy) 


O'. 


bg\xi,yx2) fg{xi,y^)fg{x2, 


+^bq^^\y^i^^2)fg{xi,y^)fg{x2,y^) + {b^b) 


+ 


47r^ 


dyidy2Pqg{yi)Pqg{y2)fg{xi,y^)fg{x2,y^)dfj^’^^\yiXi,y2X2) ) , (36) 


which is our main result. Note that in the general case in which ^ yp, the expansion 

_ 2 

Eq. (l30l) should be viewed as an expansion in powers of asiyu); the log is L = In^^; all 

'IP^b 

PDF should be evaluated at y = yp, and all hve-flavor scheme partonic cross-sections should 
be evaluated at the appropriate scale y?p). Strictly speaking, in this case 

the argument of the strong coupling in the term in Eq. (1361) which is linear in L should be 
as{yR)as{yp) = {as{yR)f{l + 0(al)). 

It is easy to see explicitly that, if the 6-PDF is expressed in terms of its values at y^ = ml 
using Eq. ([23]), the FONLL-A expression differs from the four-flavor scheme result by terms of 
order a^, namely, the difference term 

^(d) ^ ^(5),(0) ^ «^(^2)^(5),(1) ^ ^ 37 ) 

is 0{al). Indeed, Eq. (12^ implies that all contributions to the logarithmic ones are 

0{al). We then have 




dxidx2' 


fb{xi, y^)fb{x2, y^)d^’’^ ’ {xi,X2) 




dyidy2Pqg{yi)Pqg{y2)fg{xi, y^)fg{x2, {yiXi,y2X2 


(38) 


<XsdfJ’^^\xi,X2) {fg{xi,y‘^)fb{x2,y'^) + {xi X2)) 


271 


dyPqgiy) (dlf^^\xi,yx2) + alf^^\yxi,X2)) fg{xi, y^)fg{x2, y"^)] > + 0 {al). 


Substituting Eq. (12^ in Eq. (I38|) all terms in Eq. (138|) cancel, as expected 


9 











We can now stndy the phenomenological implications of onr resnlts. Leading-order fonr- 
flavor scheme predictions have been obtained nsing a modihed version of the SHERPA Monte 
Carlo generator [ 21 ] which we tested against results obtained in Ref. [ 12 ] and Ref. [ 22 ]; for 
NLO results (which we will also show for comparison) this has been further interfaced to 
the OpenLoops code [ 23 ]. Four-flavor scheme results are obtained using nj = A NNPDF3.0 
LO PDFs [ 23 ] with a^g{mz) = 0.118. Five-flavor scheme predictions are obtained using the 
bbh@nnlo code [ 13 ] with the Uf = 5 NNLO NNPDF3.0 parton set [ 23 ]- For FONLL-A, results 
for the central scale choice have been obtained in two different ways. First, we have recomputed 
the four-flavor scheme result, but now using Uf = 5 NNLO NNPDF3.0 PDFs, and we have 
combined this with our implementation of Eq. (I36|) . Then, we have checked that we get the 
same answer by combining this four-flavor scheme result with the hve-flavor scheme one from 
the bbhSnnlo code, and adding an implementation of the subtraction term Eq. ffT3ll . Scale 
variation plots have then been produced using this second combination. In all cases, the strong 
coupling provided with the PDF set has been used through the LHAPDF interface [ 23 ] • The 
b mass in FONLL expressions has been identihed with the pole mass, for which we have taken 
the value mb = 4.72 GeV; this corresponds to the MS value Wibirnb) = 4.21 GeV through the 
two-loop relation of Ref. [ 23 ], which we implemented in order to evaluate the bottom Yukawa 
coupling in the MS scheme at p = Like as and the PDFs, Yukawa couplings are evolved 
at NNLO in the hve-flavor scheme in all contributions to the FONLL expression. 

In Fig. [3]we compare the cross-section computed in the four-havor, hve-havor and FONLL-A 
scheme. Results are shown as a function of the Higgs mass. Here and henceforth, uncertainty 
bands are obtained by varying the renormalization and factorization scales fiR and fip inde¬ 
pendently by a factor of 2 about the central value iip = = discarding the two extreme 

points fiR = Afip and pf = 4//^, and taking the envelope of results. In the same hgure we also 
show the curve obtained using the so-called Santander matching of Ref. which is given by 


O'S-M 


_|_ yj (tFO 
1 + w 


(39) 


with w = Inmn/mb — 2: this reproduces the hve-havor scheme result when w —)■ (X), and the 
four-havor scheme one when to = 0. This prescription was suggested in Ref. to be used 
with the highest-order available four- and hve-havor scheme results. Here, we show it using the 
LO four-havor scheme result in order to provide a meaningful assessment of the diherences in 
comparison to FONLL-A. 

The four-havor scheme result is rather smaller than the hve-havor scheme one, and it is 
ahected by a signihcantly larger scale uncertainty, as one expects of a LO computation. The 
FONLL and hve-havor scheme results are very close, with, for mh = 125.09 GeV, the FONLL 
prediction just below the hve-havor one, with a somewhat larger uncertainty. Note that the 
four-havor scheme result shown in the plot is determined using LO PDFs, while the four-havor 
scheme result that enters the FONLL combination is consistently computed with NNLO PDFs, 
as discussed above. We have verihed that the latter would be yet lower, further away from the 
hve-havor scheme results, as one expects due to the fact the the LO gluon is typically larger. 
This shows that mass ehects for this process are small, though not negligible in comparison 
to the scale uncertainty on the hve-havor result, as we will see shortly. The fact that mass- 
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Figure 3: The total inclusive cross-section computed in the four-flavor scheme at LO (red), in 
the five-flavor scheme at NNLO (blue), and in the FONLL-A scheme (green). The Santander 
matching Eq. [39] of the four and five-flavor scheme results is also shown (purple). Both the 
absolute result (top) and the ration to the FONLL-A prediction (bottom) are shown. 


corrections at leading order are small was already noticed in Ref. Such a quantitative 

conclusion cannot be arrived at using the Santander-matched result, which simply interpolates 
between the four- and five-flavor scheme results. 

The scale dependence of the various results of Fig. [3]is shown in Fig. IHfor mu = 125.09 GeV. 
The four- and five-flavor scheme results display a significant renormalization scale dependence. 
The four-flavor scheme result drops significantly as the scale is increased because of the re¬ 
duction in value of as, while the five-flavor scheme results grows because the residual, weaker 
0{ag) dependence has the opposite sign (NNLO corrections are negative) combines with the 
growth of the Yukawa coupling with scale. Interestingly, this scale dependence cancels to a 
large extent both in the FONLL-A and Santander matched results. As a consequence, the 
mass-correct ions included in the FONLL-A result, and the scale dependence of the five-flavor 
scheme computation are of comparable size, with the FONLL result below the massless one at 
the upper range of the scale variation, and above it for lower scale choices, and specifically if 
the renormalization scale is fixed at /ir = ^s recommended in Refs. [8l[T6l[28], with a 

crossing point just below fiR = mn- 
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The factorization scale dependence is very mild in all schemes, except for FONLL, where it 
turns out that the scale dependence is of the same order as the mass-corrections, which as we 
have seen are small but not negligible. In fact, the factorization scheme dependence shown in 
the plot has been determined using as argument of the strong coupling for the term in Eq. fl36|) 
which is linear in L as{fiR)a{fiF), as discussed above. If one makes the choice which 

is equivalent up to subleading term, the scale dependence changes (and in fact it becomes 
stronger) by an amount which is comparable to the scale variation itself. This means that 
corrections of relative order ln(/ij^//ij7) to the mass-corrections are not negligible on 

the scale of the mass-corrections themselves. They could only be accounted for by upgrading 
the four-flavor scheme computation to NLO. 

Finally, in Table[I]we collect our results with mn = 125.09 GeV and /i = mn or /i = 

For comparison, in addition to the results shown in Figs. IBM] we also show the best available 
calculation in the four-flavor scheme (NLO) and its Santander matching to the NNLO hve-flavor 
result. 



(pb) 

(pb) 

^FONLL 

(pb) 

(pb) 

^S-M 

fi = mH 

U-’JO_o.o3 

n oo+U.L^b 
'^•^^-0.06 

u-oo_o.oi 

o.55^;;:^;( 

“TTocTorra 

0-56:-j^ 

fi = {niH + 2mb)/4 

0.61 

0.41 

0.82 

0.56 

0.42 

0.57 


Table 1: The total cross-section computed for mu = 125.09 GeV in the five-flavor scheme at 
NNLO, the four-flavor scheme at LO, and matching the two with FONLL-A, or with Santander 
matching (denoted as cr®“^). The NLO four-flavor scheme result, and its Santander matching 
to the five-flavor scheme are also shown for comparisons. Results are given for jj, = niH (top 
row) and /i = {mn + 2mb)/4 (bottom row). For /i = mn we also show the uncertainty band 
obtained from scale variation (see text). 

In summary, we have shown how to consistently match the four- and hve-flavor scheme 
computations of Higgs production in bottom-quark fusion. We have found that a fully matched 
computation allows detailed quantitative comparisons between the computations in various 
schemes, unlike other more phenomenological approaches. However, for competitive precision 
phenomenology, the results presented in this paper should be upgraded to include the four-flavor 
scheme result up to NLO: indeed, the factorization scheme dependence of the mass corrections 
turns out to be comparable to their size. Such an upgrade is possible by using the scheme 
presented here, in its FONLL-B version, which requires an in principle straightforward, though 
in practice somewhat laborious extension of the techniques presented in this paper: this is the 
object of ongoing work. 
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